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Tagging photon interactions at the LHC 

X. Rouby a * 

a Universite catholique de Louvain, Center for Particle Physics and Phenomenology (CP3), 
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Photon interactions at the LHC result in striking final states with much lower hadronic activity in the central 
detectors than for pp interactions. In addition, the elastic exchange of a photon leads to a proton scattered at 
almost zero-degree angle. Tagging photon interactions relies on either the use of large rapidity gaps or on the 
detection of the scattered proton using very forward detectors. The studies related to such detectors are presented, 
including their characterization, their acceptance and reconstruction performance. Limitations due to the LHC 
beamline misalignment and possible solutions are also given. 



1. Introduction 

A significant fraction oipp collisions at the LHC 
will involve quasi-real photon interactions occur- 
ring at energies well beyond the electroweak en- 
ergy scale [T|. The LHC can therefore be con- 
sidered to some extend as a high-energy photon- 
proton or photon-photon collider. The initial 
comprehensive studies of high energy photon in- 
teractions at the LHC were recently reported [21 
I3I4I5I6I7I5] . Clear identification of the remark- 
able photon-photon and photon-proton collisions 
out of the large sample of pp events requires good 
experimental tagging techniques. Two signatures 
characterizes the photon interactions: the pres- 
ence of a very forward scattered proton and of a 
large pseudorapidity region of the detector, de- 
void of any hadronic activity, usually called large 
rapidity gap (lrg), in the forward directions. 

During the phase of low luminosity (i.e. sig- 
nificantly lower than 10 33 cm~ 2 s -1 ), the event 
pile-up is negligible. Thanks to the colour flow 
in pp interactions between the proton remnant 
and the hard hadronic final states, a simple way 
to suppress generic pp interactions is to require 
LRGs. The lrg condition can be applied using a 
cut based on the energy measured in the forward 
detector containing the minimum forward activ- 
ity (3 < \rj\ < 5). More details are given in [3]. 
Tagging photon interactions with LRG could rely 
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only the central detector. But as the forward pro- 
ton is not measured, the event final state is less 
constrained. At higher luminosity, the lrg tech- 
nique cannot be used because of large event pile- 
up. Therefore the use of dedicated very forward 
detectors (vfds) is mandatory in order to retain 
pp backgrounds low. 

Both CMS and ATLAS experiments have ex- 
tended their coverage in pseudorapidities to the 
forward regions [9]. This includes near-beam 
calorimeters and tracking detectors. The trans- 
port of particles through the LHC beamlines from 
the interaction point (ip) to the forward detectors 
is simulated with Hector [10]. The simulator is 
based on a linear approach of the beamline optics, 
implementing transport matrices from the optical 
element magnetic effective length, and with cor- 
rection factors on magnetic strength for particles 
with non nominal energy. Hector deals with the 
computation of the position and angle of beam 
particles, and the limiting aperture of the optical 
elements. The measurement of the position and 
the angle of particles in dedicated near-beam very 
forward detectors, hundreds of meters away from 
the IP, helps in reconstructing the kinematics of 
the event at the related central detector. Hector 
links the information from very forward detectors 
(vfds) to the one from the central detector, by 
precise calculation of the particle trajectory. The 
measurement in the central detector of exclusive 
dimuon final states, from 77 collisions or diffrac- 
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tive photoproduction of T mesons, can be used 
for a very precise calibration of the VFDs [5] . 

2. Simulation techniques 

2.1. Physical description 

The simulator relies on a linear approach to 
single particle propagation. A beamline con- 
sists of a set of optical elements, amongst 
dipoles, quadrupoles, drifts, collimators, kickers 
and vfds. The optical elements are described by 
their magnetic field, their length and their aper- 
ture. In turn, a set of particles, with all possible 
smearings of initial positions, angles or energies, 
is propagated by Hector through the beamline, 
particle by particle. The first terms of the Taylor 
expansion of the magnetic field can be interpreted 
as dipolar (i = ^B y ), quadrupolar (k = 
and sextupolar fields. 

In the co-moving coordinate system, neglecting 
small deviations (x -C R, y <C R) and small mo- 
mentum loss (Ap <C p), this leads to the following 
equation of motions for a particle traveling along 
path length s through a magnetic clement [?] : 

\ y"(s) + k(s)y( S )=0. 

The solution (x(s),x'(s),y(s),y'(s)) to these 
equations can be expressed as a linear combi- 
nation of the initial values (xo,x ,yo,y'o), where 
the rotation matrices are defined by the proper- 
ties of the optical element (length and magnetic 
field). Each beam particle is represented by a 
phase space vector and each optical element by 
a transfer matrix by which the vector is multi- 
plied. The propagation of a single particle is thus 
the rotation of the phase space vector by the n 
optical element matrices. 



X(s) 



X(0) MiM 2 . 
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The whole beamline is modeled as a single 
transport matrix acting on each particle phase 
space vector (no intrabeam interactions). The op- 
tical element description also refers to its physical 
aperture. When a particle is propagated through 
an optical element, the compatibility between its 



trajectory and the optical element aperture is 
checked. 

2.2. Implementation 

Beam particles are described by a 6- 
components phase space vector X = 
(x, x', y, y 1 , E, 1), where (x, x') and (y, y') are the 
horizontal and vertical coordinates and angles; 
E is the particle energy. The sixth component 
is just a factor used to add an angular kick on 
the particle momentum direction. The optical 
elements (dipoles, quadrupoles, drifts, . . .) are 
modelled by 6 x 6 transport matrices: 



M, 



where 

• A (and B) blocks refer to the action (fo- 
cusing, defocusing, drift) on horizontal (and 
vertical, resp.) coordinate and angle. 

• D terms reflect the dispersion effects of the 
horizontal dipoles on off-momentum parti- 
cles. 

• K factors are the angular action of kickers. 

The chromaticity, or the energy dependence of 
the transport matrix, is implemented by rescaling 
the magnetic field terms (R, k, K) with a factor 
( )■ The propagation of particles different 
from protons is also possible by rescaling these 
magnetic field terms: 
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h = R, k, K; (2) 



where q p is the particle charge. Forward parti- 
cles from the final state can be then propagated 
through the beamline via Hector. 

3. Beamline simulation 

Knowing the optics tables for both lhc beams, 
their trajectories can be compared simultane- 
ously, in both top and side views, for the two 



Proceedings of workshop on "High energy photon collisions at the lhc' 



LHC beams aside (incoming beam 2 and outgoing 
beam 1, at the IP5 and ipI: Fig. [T]). The top 
view shows the beams on the horizontal plane, 
clearly depicting the crossing angle at the IP5, 
and the beam separation after 70 m away from 
the interaction point. The bending of the sector 
dipoles has been switched off in order to makes 
graphics more clear - this is why both beams are 
straight and parallel after 250 m. However, the 
optical elements have been shifted (without tilt) 
in the horizontal plane by the half of the beam 
separation distance, from 180 m away from the 
IP, in order to match the ideal beam path: a pro- 
ton with nominal energy and on the ideal orbit 
should travel through optical elements in their 
geometrical center. The side view in turn shows 
the beams in the vertical plane. In addition, the 
major optical elements have been drawn: rectan- 
gular dipoles in red, sector dipoles in light green, 
and quadrupoles in yellow and blue. 




s[m] 



Figure 1. Top and side views of both LHC beams 
around the ip5. The IP is located at s = m. 
Beam 2 (red) flows from the right to the left and is 
seen before its passage at the IP. Beam 1 (black) 
is seen downstream, passing from the left to the 
right after the crossing. In this graph, the bend- 
ing effect of the sector magnets (S-Dipoles) has 
been switched off, thus rectifying the beam path 
to a straight line after 250 m 10J . 



The simulation of proton beams is needed for 
the validation of the simulator results, by a direct 
comparison to the results of other simulators. 

3.1. Calibration 

The calibration of the reconstructed variables 
at forward detectors can be well maintained using 
the physics processes and the central detectors. 
At hera, for example, the elastic p meson photo- 
production was used where the momentum of the 
scattered proton could be deduced from two de- 
cay charged pions using the central tracking. At 
the LHC, the two-photon exclusive production of 
dimuon pairs and the diffractive T photoproduc- 
tion seems good calibration processes [718] . with 
large statistics for the detectors at 420 m. The 
visible cross section of exclusive dimuons is large 
(er < 7 pb), including the acceptance of central 
detectors. The measurement of the final state 
muon pair provides an estimate of the proton en- 
ergy loss x — ^/-E'beam' wn i cn can be matched 
with the proton position in the vfds. From the 
dimuon longitudinal momentum P z and invariant 
mass M, 



[ifii one gets 



V(M W )2 + (P z 



This should allow for a run-by-run calibration of 
the scattered proton energy scale within a full 
acceptance range. Finally, the expected recon- 
struction power of central detectors is excellent 
for such dimuon events, giving the proton energy 
uncertainty about 10~ 4 per event [11] . Recent re- 
sults with the full similation of the detector show 
even better resolutions, at the order 5 x 10~ 6 (de- 
pending on x) [5]. One should note, however, 
that using this process it is not possible to check 
the angular reconstruction, and that it has much 
more limited statistics within acceptance of the 
detectors at 220 m. 

4. Very Forward Detectors 

In the following, the use of very forward de- 
tectors (vfds) located at 220 m and 420 m from 
IP5 are discussed, as taggers for photon interac- 
tions. The assumed location of the first detector 
is (s = 220-224 m, x = 2000 ^m), and of the sec- 
ond one at (s = 420 - 428 m, x = 4000 ^m). No 
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hypothesis is made on their detection efficiency or 
their resolution, unless quoted. We consider here 
vfds providing 2D-measurement (x and y coor- 
dinates), each consisting in fact of two stations 
separated by 4 and 8 m as a lever arm for the 
angle measurement, with no magnetic element in 
between. 

4.1. Acceptance 

Using Hector's aperture description, it is pos- 
sible to identify the characteristics of the protons 
that will hit the vfd. The exchange of a photon, 
leaving the proton intact, results in a proton en- 
ergy loss (Eioss) and a scattering angle, directly 
linked to the four-momentum transfer squared 
(t), or equivalently to the photon virtuality (Q 2 ). 
The acceptance windows of the VFDs can be com- 
puted by performing scans in (Ei oss , t) and com- 
puting the probabilities of reaching the detectors. 
The figure [2] shows the contour plots of the detec- 
tors acceptance, in this (Ei oss ,t) plane. The VFD 
acceptances mostly depend on Ei OSSl and have a 
very small sensitivity in t, within a large t range. 
Corresponding profile at fixed virtuality is shown 
in figure [3 

The total diffractive cross-section at the lhc 
is very large, resulting in a high rate of diffrac- 
tive protons hitting the VFDs. As a result, it 
causes extremely high irradiation levels. About 
10 14 hits/cm 2 are expected per year at low lumi- 
nosity. This illustrates the need for very radiation 
hard detectors. 

4.2. Chromaticity grids 

Once the acceptance windows of very forward 
tracking detectors are defined, it is interesting to 
see matching between the proton variables at the 
IP and those measured by vfds. Depending on 
their energy and angle at the IP, forward protons 
will hit the VFDs at various positions. Drawing 
iso-energy and iso-angle curves for a set of sam- 
ple protons produces a grid in the measurement 
related variables, {xx,X2) or (x\,6). Due to op- 
tics of the LHC beamlines, the grid unfolds itself 
in a much clearer way in the latter plane, and is 
almost invisible in the former one. The energy de- 
pendence of the transfer matrices implies a defor- 
mation of the grid - without such a dependence, 
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Figure 2. vfd proton acceptance for the LHC 
beam 1 around the ip5. The vfd is located at 
(s = 420 m, x = 4000 jam). This map shows 
contours of 25%, 50%, 75% and (plain curve) 
100% acceptance. The acceptance is roughly rect- 
angular, i.e. independent of t. The angular 
kick coming from the large momentum transfers 
leads to an increasing smearing of the graph lower 
edge [ID]. 



| Acceptance at 420m (4000 urn) tor beam 1 




E,„„(GeV) 



Figure 3. vfd acceptance as a function of energy 
loss, for two fixed virtualities (Reminder: Q 2 = 
— t). See previous figure for more details [ID] . 



the grid would be a parallelogram. One should 
note, that uncertainty of the transverse position 
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of the proton vertex at the IP results effectively in 
smearing the chromaticity grids. Anyway, these 
chromaticity grids provide a straightforward tool 
for unfolding the energy and angle at the IP of the 
measured particle. The grid in figure [4] is calcu- 
lated in the energy range accessible to the VFDs. 



| Chromaticity grid at 220 m 



-200 

















: (ioo,o) 












) 


















































{Energy, 0) 














— iso-energy (i^evj 
--- Iso-angle ((irad) 




-j r (1000,500) 








.... i .... i .... i .... j. .. . 






5 1 


1 


5 


20 25 





x 1 (mm) 



| Chromaticity grid at 420 m 
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Figure 4. Chromaticity grids: iso-energy and 
iso-angle lines for the vfds at 220 m (above) 
and 420 m (below) away from the ip5, for the 
LHC beam 1. The graphs show the positions 
x\ and angle 9 X \ of protons, given the en- 
ergy loss [0; 1000] GeV and the angular kick 
of [0; 500] /irad. The energy dependence of 
the transfer matrices induces a deformation of 
the grid, worsening the reconstruction power at 
higher angles [TP"] . 



4.3. Reconstruction 

The reconstruction of the event kinematics is 
possible with Hector, from the vfd measure- 
ments. If a beam particle has exchanged a photon 
at the IP, one could reconstruct photon's energy 
(E) and virtuality (Q 2 )- The particle energy at 
a given position in the beamline is obtained from 
the measured particle position and angle within 
the matrix formalism by solving these equations: 

| i s = a s x + b s x' a + d s E 
\ x' s = a s x + s x' o + J S E 

The transfer matrix of the beamline yields the 
coefficients a, b, etc. The introduction of an en- 
ergy dependence on the strength of optical el- 
ements refines the transfer matrix, becoming a 
function of E: a s (E), b s (E),... This dependence 
will introduce non linearities. 

Neglecting the terms a s xo and b s x' leads to sig- 
nificant sensitivity to the non-nominal values of 
the average vertex position and beam direction 
(tilt) at the IP. While the average vertex posi- 
tion can be very well measured using the central 
detectors, the beam tilt is more difficult to control 
- it should be known to better than 10 — 20 /irad 
to avoid causing a significant bias. One can re- 
construct the energy and scattering angles using 
position measurements at two detector stations at 
the same time. This requires to solve the equa- 
tions for x s at both detectors for E. The a s , b s 
and d s coefficients of the transfer matrix depend 
on energy with rather complicated shapes, as they 
are products of many magnets matrices. One ef- 
ficient way to get those is then simply to fit each 
coefficient as a function of energy. Various fitting 
functions were tried, but a quadratic fit proved to 
be sufficient to avoid any visible bias or resolution 
degradation. 

Using these fitted coefficients, one can easily 
get a formula for xo as a function of x at both 
detectors and of the energy. For each pair of de- 
tector x coordinates, the used method is to nu- 
merically find the root (xo = 0) of the formula 
to get the energy corresponding to a xq — 0, and 
thus neglect the interaction point transverse ex- 
tension. 

This method allows to get reconstructed energy 
independently of the angle of the particle at the 
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IP. Energy reconstruction resolutions does not 
degrade anymore if the particle transverse mo- 
mentum pt rises, but the price to pay is that 
the detectors resolutions become critical. In par- 
ticular, it means that the uncertainty of the re- 
constructed angle at a given detector location 
should be better than the beam angular diver- 
gences there (6 and 1.5 /zrad, respectively). The 
transverse momentum can thus also be computed 
from the matrix coefficients once the energy has 
been reconstructed. 

The reconstructed energy resolution stays very 
good even with non- negligible initial particle px- 
Figures [5] and [6] show the px and energy depen- 
dence of resolutions on the transverse momentum 
and energy, for various detectors resolutions. As 
expected, the energy resolution is independent of 
Pt but is sometimes sensitive to the energy as 
expected from the chromaticity grids of Fig. |4] 
The effect of the detector resolution, which was 
absent for the simple reconstruction method de- 
scribed previously, is now clearly visible even for 
an excellent resolution of 5 /im, especially for de- 
tectors at 220 m from IP. 

In summary, the proton angular distribution 
at the IP affects the most the energy reconstruc- 
tion. In contrast, the vertex lateral distribution 
has negligible impact. As a result, the most im- 
portant beam parameter, which can change run- 
to-run, is the beam tilt at the IP. In principle, 
it can be indirectly controlled by the beam po- 
sition monitors (bpms) at 220 and 420 m, but 
more direct tilt measurements are favored, as by 
using BPMs next to the IP, or by monitoring the 
direction of neutral particle (photons, neutrons, 
or neutral pions) production in the zero-degree 
calorimeters (zdcs). 

The energy resolution squared of the scattered 
proton (Te can be then approximately decom- 
posed into four terms: 

a \ = ^0 + °ltx + a l.ng + °def, 

where the nominal beam energy dispersion a® ~ 
0.8 GeV, the contribution due to the vertex 
spread a vtx ~ 0.7 GeV at 420 m and 1-2 GeV 
at 220 m, the contribution due to the detector 
resolution, neglecting the angular effects, a^ et is 
small for resolutions better than 50 /im, and the 



contribution due to the proton non-zero angle at 
the IP, a ang , which is very sensitive the angular 
reconstruction at the VFDs; even for a zero-degree 
scattering if one neglects effects due to the beam 
angular divergence at the IP, a ang w 1 GeV at 
420 m and w 3 GeV at 220 m. 
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Figure 5. Resolution of the reconstruction of the 
particle transverse momentum pr, as a function 
of the energy loss (above) and of the transverse 
momentum (below), for vfds at 420 m from the 
IP5. Dots correspond to different scenarios of de- 
tectors resolutions, namely perfect detectors (cir- 
cles), 5 /im (squares) and 30 /im (triangles) spa- 
tial resolution [10] . 



5. Misalignments 

The misalignment of lhc optical elements 
could have a significant impact on the measure- 
ments with very forward detectors. As the de- 
flection of the particle paths depends on their po- 
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Figure 6. Resolution of the reconstruction of the 
particle energy loss E, as a function of the en- 
ergy loss (above) and of the transverse momentum 
(below), for VFDs at 420 m from the ip5. Dots 
correspond to different scenarios of detector reso- 
lutions, namely perfect detectors (circles), 5 /xm 
(squares) and 30 /im (triangles) spatial resolu- 
tion [TU]. 



sitions in quadrupoles, a misplacement of these 
optical elements implies a change in the nominal 
beam position. In turn, as the accurate position 
measurement with the forward tracking detectors 
(as well as the information inferred from the seg- 
mentation of forward calorimeters) is referred to 
the ideal beam location, changing this reference 
results in a biased reconstruction of the measured 
particles. 

Figure [7] shows the impact of possible shifts 
(0.5 mm) of the beamline quadrupoles on the en- 
ergy reconstruction with vfds at 420 m. The 
reconstruction assumes ideal beamline in which 
only one quadrupole at a time is separately 
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Figure 7. The misalignment of the LHC 
quadrupoles bias the energy reconstruction. The 
graphs show the bias for the reconstruction (with 
the trivial method) of a 100 GeV energy loss 
assuming the misaligned quadrupoles at various 
positions. Each element is separately shifted 
(500 /xm) , assuming a perfect alignment for the 
rest of the beamline. The impact of the misalign- 
ment can be important. Even a perfect knowledge 
of the actual beam position at the vfd (below) 
does not compensate for this bias, depending on 
the position of the misplaced quadrupole [TO] . 



moved. Effects even higher than 10 % could be 
expected. Similar results are observed due to the 
beam tilts at the IP [10] . One can partially cor- 
rect for these effects using information from the 
beam position monitors, but better results are ob- 
tained using a physics calibration process like the 
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two-photon muon pair exclusive production or the 
diffractive photoproduction of T, at least in case 
of the vfds at 420 m [5]. 

These misalignment effects and the corrections 
are illustrated (Fig. [8]) by the study of two- 
photon exclusive production of the SM Higgs bo- 
son (pp -> pp(77)H) with M H = 115 GeV. The 
measurement of the energy of two scattered pro- 
tons straightforwardly yields the boson mass, by 
means of the so-called missing mass method. As 
a consequence, an uncorrected measurement with 
misalignment leads to a bad mass calculation. 
A quadrupole (mqxa1r5, s — 29 m) close to 
the IP is shifted by 500 /im. The misalignment- 
induced change in the vfd acceptance is visible. 
The limitations of the beam-position-based cor- 
rections are clearly visible, even assuming no sys- 
tematic errors, while the muon-calibration stays 
unbiased (though only a relatively small sample 
of 700 dimuon events was used to get the correc- 
tion factors). 

6. Summary and outlook 

Photon interactions at the LHC result in strik- 
ing final states with much lower hadronic activ- 
ity in the central detectors than for pp interac- 
tions. While forward rapidity gaps can be used 
for their tagging in a low pile-up environment, the 
use of very forward detectors is mandatory for lu- 
minosity starting at L — 2 x 10 33 cm _2 s -1 . The 
good simulation of particle transport in the LHC 
beamline is very important. It allows to charac- 
terize the detectors and to reconstruct the event 
from the scattered proton measurement. Impact 
of beamline misalignement has been studied as 
well as the calibration possibilities from the ob- 
servation of exclusive dimuon final states. 
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